*******  replication Table A13: logit models

eststo clear


** MNLF
use "MNLF_Target_Final.dta" , clear
    


quietly: logit event ramadan christian_4days christian_2days christian_1days nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day  rain hot lag_event cubic1-cubic3 i.region  , vce(cluster region) from (ramadan=-0.9)

eststo model1


** MILF: clogit does not work; ombinations results in numeric overflow; computations cannot proceed

use "MILF_Target_Final.dta" , clear


quietly: logit event ramadan christian_4days christian_2days christian_1days nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day  rain hot lag_event lag_gov_counter cubic1-cubic3 i.region  ,  vce(cluster region) from (ramadan=-0.5)
eststo model2



** BIFF
use "BIFF_Target_Final.dta" , clear



quietly: logit event ramadan christian_4days christian_2days christian_1days nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day  rain hot lag_event lag_gov_counter cubic1-cubic3  i.region  , vce(cluster region) 
eststo model3



*** Patani before 2000
use "Patani_bef_Target_Final.dta" , clear




quietly: logit event ramadan  Bud_holiday nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day rain hot  lag_event cubic1-cubic3  lag_gov_counter  i.region,    vce(cluster region) from (ramadan=-0.09)
eststo model4


*** Patani after 2000
use "Patani_aft_Target_Final.dta" , clear




quietly: logit event ramadan  Bud_holiday nonreli_holiday  /*
*/ nation_election Islamic_3day  Islamic_2day Islamic_1day rain hot  lag_event cubic1-cubic3 lag_gov_counter   i.region, vce(cluster region) 

eststo model5



*** GAM
use "GAM_Target_Final.dta" , clear


quietly: logit event ramadan christian_4days christian_2days christian_1days nonreli_holiday  /*
*/ nation_election Islamic_3day  Islamic_2day Islamic_1day  rain hot lag_gov_counter lag_event  cubic1-cubic3 i.provincode, vce(cluster provincode)  from (ramadan=-0.1)
eststo model6


esttab using Table_A13.tex,  order(ramadan Islamic_3day Islamic_2day Islamic_1day  christian_4days christian_2days christian_1days Bud_holiday    nonreli_holiday  nation_election lag_gov_counter  rain hot attack_gov_region_lag cubic1 cubic2 cubic3 ) label star(+ 0.1 * 0.05 ** 0.01 *** 0.001) se b(3) replace booktabs ///
  alignment(D{.}{.}{-1})      mtitles("Model 1" "Model 2"  "Model 3" "Model 4"  "Model 5" "Model 6")                   ///
  title(Fixed Effects Model table\label{tab1})

